clear all;
clc;

filename = 'sub_bp_dr_phi_100axons_167nodes_9um.txt';
drphi = load(filename);
[m,n] = size(drphi);

drphi(26,:) = drphi(27,:);
drphi(43,:) = drphi(44,:);

ck = (drphi==drphi);
[ri,ci] = find(ck==0);

ri
ci

dlmwrite(filename,drphi,'delimiter',' ');
